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The production of both charged and neutral electroweak gauge bosons in association 
with b jets has attracted a lot of experimental and theoretical attention in recent years 
because of its central role in the physics programs of both the Fermilab Tevatron and 
the CERN Large Hadron Collider. The improved level of accuracy achieved both in the 
theoretical predictions and experimental measurements of these processes can promote 
crucial developments in modeling 6-quark jets and 6-quark parton distribution functions, 
and can provide a more accurate description of some of the most important backgrounds 
to the measurement of Higgs-boson couplings and several new physics searches. In this 
paper we review the status of theoretical predictions for cross sections and kinematic 
distributions of processes in which an electroweak gauge boson is produced in association 
with up to two 6 jets in hadronic collisions, namely pp,pp V lb jet and pp,pp 
V 2b jets with V = , Zf'y* ,'y. Available experimental measurements at both the 

Fermilab Tevatron and the CERN Large Hadron Collider are also reviewed and their 
comparison with theoretical predictions is discussed. 

Keywords: Quantum chromodynamics; electroweak gauge bosons; bottom quarks; 
proton-proton inclusive interactions 

PACS Nos.: 13.85.Ni, 12.38.Bx, 14.65.Fy, 14.80.Bn 


1. Introduction 

The production of both charged (W^) and neutral {Z/j*,j) electroweak gauge 
bosons in association with b jets (generically denoted as V + b-jet production) has 
attracted a lot of experimental and theoretical attention in recent years because 
of its central role in the physics programs of both the Tevatron and the Large 
Hadron Collider (LHC). Many formal developments have improved the accuracy 
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of theoretical predictions to the level of accuracy necessary to compare with the 
results that have meanwhile been presented by both the Tevatron and the LHC 
experiments. In this paper we present an updated review of this topic and collect 
the most important results that can serve as reference for future studies. 

From a physics standpoint W and Z production with 6 jets is first of all one of the 
main irreducible backgrounds to the measurement of a Higgs boson in the WH and 
ZH associated-production channels, which played a leading role for Higgs-boson 
searches at the Tevatron and will be instrumental in constraining the couplings 
of the discovered Higgs boson^’^ at the LHC. W + b jets is also and important 
irreducible background in the measurement of single-top production, and Z -\-b jets 
can mimic the H —)■ ZZ* —)■ Al (^=lepton) final state when associated B mesons 
decay leptonically. In general, electroweak vector bosons (namely and Z) and b 
quarks often appear in the decay chain of heavy massive particles in models beyond 
the Standard Model (SM), making of the V -I- 6 -jet processes important backgrounds 
for searches of new physics at the LHC. 

At the same time, the production of a SM vector boson {W^, Zj^*, and 7 ) 
with heavy quarks (both b and c) is per se an important test of QCD for two main 
reasons. The first reason is that it allows to study the dynamics of b (and c) jets and 
improve their modeling in many other more challenging SM processes that represent 
important backgrounds to new physics searches (consider for example tt + b jets). 
The second important reason is that the production of , Z, and 7 with heavy 
quarks (b and c) offers a unique possibility to measure the bottom- and charm-quark 
parton distribution functions (PDF) using only collider data, as will become clearer 
(for the b case) from a closer inspection of the parton-level processes that are at the 
core of H -I- 6 -jet production. 



(a) 



(b) 


Fig. 1. Examples of tree-level Feynman diagrams contributing to + b signatures at hadron 
colliders, without and with contributions from initial-state b quarks. 


The parton-level structure of V -|-6 jets differs for charged (V = ) and neutral 

(V = Zj^*,^) gauge bosons. Furthermore, such structure depends on the number 
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of light flavors {Nf = A or Nf = 5) that are allowed in the initial-state hadrons 
{p and p). Given the non-negligible mass of the bottom quark (mb ^ (^qcd), the 
question of considering a massless bottom quark with an initial-state parton density 
in the (anti)proton is non trivial and is intuitively related to the energy scale of the 
process considered. Since typically V + &-jet processes involve a range of energies 
from few tens to few hundreds of GeV, the answer is not clear cut and the problem 
cannot be addressed without considering the perturbative QGD structure of the 
corresponding physical observables (total and differential cross sections), as we will 
discuss in more detail in Section 2. Considering only the first-order (or tree-level) 
processes, let us remind the reader that + 2b jets can proceed only through 
qq' —)■ W^bb (see Figure l.a), while + lb jet can proceed through the same 
channel if we assume Nf = A and also through bq —> W^b + q' (see Figure l.b) if 
we assume Nf = 5. Similarly, Z/y -|- 2b jets can only proceed (at lowest order in 
Os) through qq —> Z/jbb and gg —)• Z/jbb (see Figure 2.a), while Z/y -|- lb jet can 
proceed through the same processes if we assume Nf = 4, or through bg —)■ hZj^ if 
we assume Nf = 5 (see Figure 2.b). Notice that in this discussion we keep only the 
Cabibbo sector of the flavor-mixing CKM matrix and neglect processes initiated by 
two b (anti)quarks (like e.g. bb —> jbb). These approximations have effects at the 
1% level in the observables studied. The diagrams in Figures l.b and 2.b clearly 
illustrate the possibility of measuring the 6 -quark PDF by measuring V + 6 -jet 
processes. 

On the theoretical side, total and differential hadronic cross sections including 
Next-to-Leading Order (NLO) QGD corrections have been first calculated for W^bb 
and Zbb production in the limit of zero bottom-quark mass,^’^ followed more re¬ 
cently by fully massive calculations {mt ^ 0) for 1 F^ 66 ,®“^ Zbb,^'^ and y 66 ® produc¬ 
tion. The NLO QGD calculations of W^bb hadronic production (for mb ^ 0), includ¬ 
ing spin-correlation effects from W —>■ Ivi decay,^ and of Zl'f*bb (with mb = 0 and 
Z/y* —>■ Nl~) are available in MCFM.^° These NLO QGD corrections can now also 
be calculated with several automated tools like for example MadGraph5_aMC@NL0,^^ 
and SHERPA^^ in conjunction with GoScun^^ or DpenLoops.^^ A study of W^bb, and 
Z/'y*bb production at hadron colliders, including NLO QGD corrections, mb effects, 
and spin-correlation effects from W and Z decays, as well as the interface with 
parton-shower Monte Carlo generators like PYTHIA^^’^® and HERWIG^^’^® via the 
MCSNLO method,^®’has been presented in Ref. 21. The NLO QCD calculation of 
pp,pp —> W^bb from Refs. 5 and 6 has also been interfaced with Pythia and Herwig 
using the POWHEG method,and is available in the POWHEG BOX framework. 
Recently the NLO QCD hadronic cross section for W^bb + j has been implemented 
in the POWHEG BDX^^ (see also the calculation of the corresponding 0{as) virtual 
corrections presented in Ref. 26). 

Higher-order QCD corrections are also known for the Nf = 5 processes. In 
particular the first complete calculation of hadronic production with at least one 
6 jet {pp,pp —>■ IF-I- 6 -I-A), including NLO QCD corrections and full 6 -mass effects. 
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Fig. 2. Examples of tree-level Feynman diagrams contributing to Z /7 b signatures at hadron 
colliders, without and with contributions from initial-state b quarks. 


has been presented in Ref. 27, updated in Ref. 28, and is now available through 
MCFM.^° The calculation was based on separate contribntions from Refs. 29, 5, 6 , and 
7. NLO QCD corrections to Z hadronic production with one b jet (pp,pp —^Z + b) 
and with at least one 6 jet {pp^pp —)■ Z + b+X) have been calculated in Refs. 30 and 
31 respectively, and are available throngh MCFM. Finally, the hadronic production 
of a photon and one b jet (pp,pp —>■ 7 + &) at NLO QCD has been calculated in 
Ref. 32 and independently cross-checked in Ref. 9. Some of these processes are now 
being computed using automated NLO QCD parton-shower Monte Carlos like e.g. 
MadGraph5_aFIC@NL0,^ although some aspects of the parton-shower generation in the 
Nf = 5 scheme still need to be studied, as we will discuss in Section 2. 

On the experimental side, both the CDF and DO experiments at the Tevatron 
as well as the ATLAS and CMS experiments at the LHC have provided a variety 
of results for V + 16-jet and V + 26-jet production, both with and without extra 
light jets. At the Tevatron, the cross section for + 16 jet has been measured by 
both CDF^^ and CDF has also measured the cross section for Z +lb jet,^® 

and DO the ratio of the cross sections ior Z + b jets and Z + j (j=light jets). The 
prodnction of^^ 7 -F 16 jet has further been measured by CDF^® and D 0 ,^®’^° as 


^ Results from MadGraph5_aMC@NL0 have been published in ATLAS and CMS papers as we will 
summarize in Section 3 
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well as of 7 + 2b jets by D0.^° At the LHC, we now have measurements of + lb 
jet and + 2b jets from ATLAS, _|_ 2 ^ from CMS,^^ as well 
as of Z + 16 jet and Z + 25 jets from both ATLAS^'^AS grid We will 

summarize the most relevant aspects of these experimental studies including the 
exact definition of the signal signatures in Section 3. 

The rest of the paper is organized as follows. We briefly review in Section 2 
the main structure of the V + 5-jet cross sections, including the first order of QCD 
corrections, for both Nf = 4 and Nf = 5. In Section 2 we also discuss the uncer¬ 
tainties associated with the theoretical predictions for these processes, and present 
more specific theoretical issues, like photon isolation, which is relevant for 7 -|- 5-jet 
production. A selection of the most recent experimental results will be presented in 
Section 3 with particular emphasis on the comparison between theory and exper¬ 
iments as presented in each of the selected measurements. Section 4 contains our 
conclusions and outlook. 

2. Theoretical Framework 

In this Section we review the most important aspects of the calculation of the 
hadronic cross section for H -|- 15 jet and 14-1-25 jet production {V = , ZjY,l) 

at NLO in QCD. In particular we focus on the meaning of presenting results in the 
four- {Nf = 4, denoted as 4F) and five-flavor {Nf = 5, denoted as 5F) schemes, on 
the matching of the corresponding parton-level calculations with existing parton- 
shower Monte Carlo generators, and on the assessment of the theoretical uncer¬ 
tainties associated with different stages of the calculation. We will substantiate our 
discussion with results taken from various theoretical studies presented in the liter¬ 
ature, and refer to the corresponding sources for further details. Several important 
aspects of existing theoretical studies are illustrated using the experimental results 
collected in Section 3. 

2.1. Four- and five-flavor schemes 

Four- and hve-flavor schemes differ by the number of partons that are allowed to 
have a non-zero parton distribution function in the colliding hadrons. Since the 
bottom quark {mb ^ ^qcd) lives in the perturbative regime of QCD, the issue of 
dehning a bottom-quark parton density is non trivial, and is intuitively related to 
the energy scale of the process considered. It is important to understand the relation 
between the two schemes and emphasize that they are just different prescriptions 
to calculate the same physical processes {V 5-jet hadronic production in our case) 
in perturbative QCD. One can indeed tailor the arrangement of the perturbative 
expansion of the cross sections to specific energy regimes, where some classes of 
corrections may become more important than others. Since 4F and 5F correspond 
to different rearrangements of the perturbative expansion of the cross section for a 
given process, the two prescriptions share more terms the higher the perturbative 
order, and in some cases significant differences can still appear at NLO in QCD. A 
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recent systematic study of 4F versus 5F schemes has been presented in Ref. 49, to 
which we refer for more details. 

V + 5-jet hadronic processes offer indeed a broad range of examples and very 
effectively lend themselves to illustrate the range of validity of fixed-order calcu¬ 
lations using either the 4F or the 5F approaches. Given the multi-scale nature of 
V + 5-jet processes, the accurate theoretical description of these processes requires 
to consider the perturbative QCD structure of the corresponding physical observ¬ 
ables (total and differential cross sections) in detail. We will return to this point 
with specific examples after having defined more in detail how 4F and 5F schemes 
are implemented in the actual calculation of R -I- 5-jet hadronic production. 

In the 4F scheme, the initial-state quarks are constrained to be the four lightest 
quarks (taken to be massless), there are no bottom quarks in the initial state, and 
the final-state bottom quarks are considered massive. Both V + 25-jet and V + 15- 
jet processes receive contributions at lowest order in QCD from qq and gg initiated 
subprocesses, specifically from qq' —> W^bb in the case of , and qq —>■ Zj^bb 
and gg ^ Z/^bb \ii the case of Z and 7 . Cross sections are computed at fixed order 
in QCD (LO, NLO, etc.) and the number of 5 jets in the final state is selected 
applying specific kinematic cuts on the transverse momentum (p^'’*** >p™*") and, 
possibly, the rapidity/pseudorapidity (|j/*'j®*| or < 7 ^™“^) of the final- 

state 5 jets, where it is understood that a given jet algorithm is in place to select 
the final-state 5 jets and light jets.*^ At NLO in QCD the cross section for R55 
production receives contributions from both the lowest order processes listed above 
and the corresponding 0{as) virtual and real corrections, more specifically: 

• one-loop virtual corrections to qq' —)■ 1F^55 and qq^gg —t Z/"fbb, respectively; 

• one additional parton real emission processes, i.e. qq' —> IF^55 -I- g and qg —>■ 
IF^55 -I- q' in the case of W^bb production; as well as qq,gg —> Zj^lh g and 
qg — > Z/”fbb -|- g in the case of Zj^bh production (notice that for convenience we 
denote by qg also the corresponding qg channels). 

In the 5F scheme, the definition of a bottom-quark PDF arises because the 
integration over the phase space of non-identified final-state bottom quarks (as in 
the case of R -I- 15-jet processes) gives origin to potentially large logarithms of the 
form: 

Ab = lnf4) - (1) 

\miJ 

where ^ is a scale related to the upper bound on the integration over the transverse 
momentum, pT^ of the final-state bottom quark,while the dependence on mb comes 
from the limit of producing a b quark with no transverse momentum. This happens 

^We notice that in a parton-level calculation, even at NLO in QCD, these jets have very limited 
structure. 

^The symbol p is often used to suggest the possibility of setting this upper bound to be the 
factorization scale of the bottom-quark PDF to be defined in the next paragraph. 
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when a final-state bottom-quark pair is produced via g ^ bb splitting of an initial- 
state gluon, and corresponds to the collinear singularity that would be present in 
the case of a gluon splitting into massless quarks (wf, —>■ 0). For scales ^ of 0{Mv) 
or larger these logarithms can be quite large. Additionally, the same logarithms 
appear at every order in the perturbative expansion of the cross section in the 
strong coupling, as, due to recursive gluon emission from internal bottom-quark 
lines, as well as virtual corrections. These logarithms could hinder the convergence 
of the perturbative expansion of the total and differential cross sections. In the 5F 
the convergence is improved by introducing a perturbatively-defined bottom-quark 
PDF,®°“^^ defined at lowest order in as as, 

b{x,^i) = J yCjs , (2) 

where g{y,n) is the gluon PDF and Pqg is the Altarelli-Parisi splitting function for 
qq given by, 

+ . (3) 

The Ah logarithms are subsequently resummed through the DGLAP evolution equa¬ 
tion,such that contributions proportional to (asAf,)” can be absorbed, to all 
orders in n, into a leading-logarithm bottom-quark PDF, while subleading loga¬ 
rithms are recursively resummed when higher-order corrections are considered in 
the DGLAP equation. With the use of a bottom-quark PDF, the 5F effectively 
reorders the perturbative expansion to be one in Ug and Ah, or agA},. In those 
kinematic regimes where a^Ah is not a small expansion parameter, the use of a 
bottom-quark PDF should improve the stability of the total and differential cross 
sections. As observed in Ref. 49, however, the argument of the Ah logarithms, as it 
arises from the collinear limit ot g —>■ bb in 4F processes, is often not a constant, 
but contains some dependence on the momentum fraction carried by the gluon {y 
in Eq. 2), which can considerably reduce the relevance of such logarithmic terms 
at each perturbative order. General arguments to estimate their relevance at all 
perturbative orders are given for 2 —> 1 and 2 —> 2 processes in Ref. 49. 

In the 5F, while the production of V -I-2& jets mainly receives contributions from 
light-quark- and gluon-initiated processes (processes like bb —>■ Zbb are negligible), 
the case of IT^-|-I6 jet and Z/q-l-lb jet have rather different characteristics. Z/y+lb 
jet starts at LO with bg —)■ Zjyh, where, according to our previous discussion, the 
initial state Aquark density has been introduced to resum the large logarithms 
arising in the calculation of gg —)■ Zjyhb, when one of the b quarks in the g ^ bb 
splitting gives origin to a 6 jet that is not tagged. At lowest order it corresponds to 
b in Eq. (2). At NLO the following 0{as) corrections are then to be added: 

• one-loop virtual corrections to —> Zjyb (notice that for convenience when we 

write h{g) in the initial state we mean generically processes initiated by both hly) 
and b{q)); 
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• one additional parton real emission processes: bg —> Zj^h + 5 , 6(7 —)■ Zj^b + q, 
qq - 7 - Zfjbb, and gg —)■ Z/^bb. 

On the other hand, in + 16-jet production the first 6 -initiated subprocess starts 
as a 2 —> 3 process {bq —>■ W^bq') because the first gluon-initiated subprocess is 
higher order {qg —> W^bb + q'). Therefore, there are two competing LO subprocess, 
namely qq' —)■ W^bb and bq —> W^bq' and the NLO QCD cross section includes 
the following 0{as) corrections: 

• one-loop virtual corrections to qq' —> W^bb] 

• one-loop virtual corrections to bq —>■ W^bq'; 

• one additional parton real emission processes: qq' —> W^bb + g, qg —> W^bb + q', 
bq W^bq' + g, and bg -> W^bqq'. 

One can notice that the subprocesses and corrections contributing to the 4F NLO 
cross section for + 16-jet production are a subset of those contributing to the 
5F NLO cross section, while this is not true for Z /7 -|- 16 jet. In the case of + 16 
we can then state that the most complete NLO QCD prediction is the 5F one, while 
in the case of Z/j + 16 jet it very much depends on the relevance of the logarithms 
that are resummed in the 5F approach as compared to the non-logarithmic terms 
that are neglected in the 5F but not in the 4F. 

We also notice that special attention needs to be paid in combining tree-level 
and higher-order corrections in order to avoid double counting of logarithms that 
appear both in a fixed-order process (like, e.g. gg —> Zbb) from gluon splitting, and 
in the corresponding 6 -initiated subprocess (like, e.g. bg —> Zb) from the expansion 
of the 6 -quark PDF. Double counting is avoided by introducing a subtraction term 
that removes from the the fixed-order cross section of a 6 -initiated process the 
logarithmic term(s) already resummed in the 6 -quark PDF. For example, the NLO 
QCD cross section for Z -|- 16 jet in the 5F should be calculated as, 

cr^F^°(Z -t 16) = cr^° -t - J dxidx2b{xi, g)g{x2, g.)cr^^ , (4) 

where the last term is the NLO subtraction term which avoids double counting of 
the Ab logarithm already included in the lowest order 6 PDF given in Eq. (2). 

It is clear by looking at the different perturbative structure of the calculation 
of IF ^66 vs. Z /766 production that we can expect 4F or 5F predictions to give 
a more accurate prediction of these processes depending on the energy (Tevatron 
vs. LHC) and kinematical regime. Apart from more specific situations on which we 
will comment in Section 3 when we discuss the comparison of existing experimental 
measurements to the corresponding theoretical predictions, we can conclude this 
synopsis of 4F vs. 5F by emphasizing some general criteria. First of all, it is worth 
emphasizing one more time that V + 26-jet processes are very well described in the 
4F at NLO in QCD. This is clear by a simple inspection of the contributing parton- 
level processes and is supported by the measurements of IF -I- 26 jets (see Table 2 
and Figure 3), Zj^* + 26 jets (see Table 4 and Figures 5 and 6 ), and 7 + 26 jets (see 
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Figure 9). On the other hand, for V + 15-jet processes 4F and 5F approaches can 
be complementary. In general, a 5F approach lends itself better to the calculation 
of more inclusive quantities (like total cross sections) and, in this case, it often 
allows to push the perturbative fixed order slightly higher since, in most cases, it 
involves simpler processes. At the same time a 4F approach can be preferable for 
the calculation of more exclusive observables, like processes with more than a single 
heavy jet in the final state. 

There are however cases in which one approach can be clearly preferable than the 
other beyond general considerations. For instance, processes that, at a given energy 
or in a given kinematical regime, are determined at LO by gg-initiated subprocesses 
will be less sensitive to large logarithms Af, that have their origin in the integration 
over the low-px region of 5 quarks produced in (/ —> 55 initial-state splitting. This 
is for instance the case of IF -I- 5 jet, where the sensitivity to A^-like logarithms 
is indeed a NLO effect (see results in Table 2 and Figures 3 and 4), or the case 
of 7 + 5-jet production at the Tevatron, where the gg-channel dominates over the 
(/g-channel production (see results in Figures 7, 8 and 9). 

Ideally, a 5F calculation with sufficiently high orders of QCD corrections included 
would give the best of both worlds, since it would include both enough fixed-order 
non-logarithmic terms (also present in the corresponding 4F calculation) and the 
resummation of several orders of potentially large initial-state logarithms. The more 
perturbative orders are included the broader is the overlap between 4F and 5F 
predictions for a given cross section. It is in general a good sanity check to assure 
that the 4F and 5F predictions confirm the existence of such overlap by showing 
more and more compatibility in going from LO to NLO, or, if needed, to higher 
order. 


2.2. Theoretical uncertainties 

Theoretical predictions are affected by both perturbative and non-perturbartive un¬ 
certainties. On one hand, hard-scattering cross sections are inherently approximate 
since they are calculated at a given order in perturbative QCD, while the algo¬ 
rithms adopted to implement the following radiation of partons (parton showering) 
are based on the approximation of accounting only for those phase-space regions 
that most contribute to extra-parton radiation. On the other hand, all hadroniza- 
tion models adopted to describe the formation of hadronic bound states after the 
last stages of parton showering contain some level of arbitrariness. 

At the parton level, the main sources of theoretical uncertainty are the depen¬ 
dence on the renormalization and factorization scales {pr, Pf) of the NLO QCD 
cross sections as well as their dependence on the choice of PDF and a^. The de¬ 
pendence on the choice of m?, can also be sizable and comparable at times to the 
systematic uncertainty from PDF and Og- This has been documented for instance 
in some IF + 5-jet studies,’^’but is otherwise usually neglected in experimental 
studies. Looking at the theoretical predictions reported in Tables 2 and 4, we can see 
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that uncertainties due to scale variation are typically at the 10-25% level, depending 
on the process, while the uncertainty from PDF and as is closer to 2-8%. On the 
other hand, from the entry in Table 2 we see that DO estimates the uncertainty due 
to mb at the 5% level. 

It is difficult to summarize the scale dependence of different V -I- 6-jet processes 
in a way that could fit the various experimental frameworks in which they are 
measured, since the behavior of the theoretical predictions depends in part on the 
specific signatures required by the experiments. There are however some general 
characteristics that only depend on the intrinsic QCD structure of these processes 
and we will review them here. In Section 3 we complement this information by 
collecting in Tables 2 and 4 existing experimental results and the quoted theoretical 
predictions, including their uncertainties, and we will add comments on the specific 
of each case if necessary. 

From several general studies appeared in the literature over the last few 
years®“®’^^’^® we can see that the cross section for all P -I- 2 6-jet processes shows 
a large residual scale dependence even after fnll NLO QCD corrections have been 
included. As explained in Section 2.1 the calculation of P -f 2 6 jets is a 4F-only 
result, and the residual scale dependence is mainly due to the opening of a new qg 
channel {qg —)■ P66 -I- q^''>) at NLO. This subprocess contributes to the 0{as) real 
corrections of P66 production and, as a tree level process, affects the cross section 
with the large scale uncertainty of a LO process. This is most visible in the case of 
IP=*=66 production, since at LO W^bb is only qq'- not gluon-induced. At the LHC, 
where the gluon density dominates over the light-quark densities, gg-initiated pro¬ 
cesses are enhanced compared to qq processes and the impact of the gg-initiated 
0{as) corrections induce a large scale dependence of the NLO cross section. This 
effect is milder, although still very visible, for Z/j + 2 6-jet production, since the 
corresponding LO processes already contain a g( 7 -initiated channel (see Figure 2), 
and are therefore less affected by the opening at NLO of a qg channel. This has 
been observed and carefully discussed in Refs. 5, 6, and 8. 

On the other hand, NLO theoretical predictions for P -|- 16-jet production show 
in general a more stable behavior with respect to renormalization- and factorization- 
scale variations when calculated in the 5F, since the logarithmic resummation in¬ 
troduced via the 6-quark PDF also resums scale-dependent logarithmic corrections 
and, as a side effect, improves the scale dependence of the cross section. On the 
other hand, the 4F predictions of P -I- 16-jet production show the same large scale 
dependence as the corresponding P -|-26-jet results. This was studied in Refs. 27, 28 
for IP + 16-jet production and can also be seen for the case of Z + 16-jet production 
in Table 4 where results for Z + lb jet production in both the 5F and 4F schemes 
are reported. 

The second major source of theoretical uncertainty at parton level is the depen¬ 
dence on the PDF with which the parton-level cross section is convoluted to produce 
the hadronic cross section. In the 4F the systematic uncertainty from PDF is due to 
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the accuracy of light-quark and gluon parton densities. In the 5F the same is true 
since the bottom PDF is defined in terms of the gluon and light-quark PDF. In both 
cases, the uncertainty induced on total and differential cross sections is of the order 
of 2-8% (see Tables 2 and 4), smaller than the uncertainty from scale dependence, 
and it will probably further improve in the future. Still, the theoretical uncertainty 
from PDF and as has not been systematically and consistently implemented in 
all experimental studies. Many results are still compared to theoretical predictions 
that use specific PDF (without comparing to others), or do not distinguish between 
differences due to the choice of Og and of the PDF set. 

In particular, the dependence on the 6-quark PDF needs to be addressed more 
carefully. Indeed we need to remember that the V -|-6-jet processes can play a leading 
role in measuring the 6-quark PDF, providing the first such direct determination 
and therefore testing the range of validity of a purely perturbative definition of the 6- 
quark PDF in terms of gluon and light-quark parton densities. It is therefore crucial 
to determine the precision with which the 6-quark PDF itself can be extracted from 

V -f 6-jet measurements. A careful determination will have to consider all available 

V -j- 6-jet processes and try to select a sample of measurements that provide enough 
information, when compared with theoretical predictions, to disentangle light-quark, 
gluon, and 6-quark initiated processes. Typically, separating samples with two 6 jets 
from samples with one 6 jet, as well as disentangling single 6-jet events in which 
the 6 jet contains two almost collinear 6 quarks or a single 6 quark can be very 
helpful. To this extent some theoretical studies^^’^® have reported predictions for 
each sample separately, and recent experimental analyses, like Ref. 56, accounts for 
the study of techniques developed to tag so called merged (66) jets. 

Finally, it is important to mention one last source of parton-level theoretical 
uncertainty, namely the effects due to double-parton scattering,®"^ that has been 
only recently formally included in the theoretical prediction for -h 26 jets and 
Zj'j* -I- 26 jets. Double-parton scattering, where aV -\-2b {V = W^,Zj"f*) parton- 
level final state is produced as the result of two scattering processes (gg^'^ —> V and 
qq, gg —> 66), can give a substantial contribution when the constituent subprocesses 
have very large cross sections (as is the case for dijet or 66 production). The factor¬ 
ization on which this estimate is based, however, is a rather crude approximation 
and relies on the measurement of the effective-cross-section parameter (ce//) for 
hard double-parton interactions. The large uncertainty that still affect such mea¬ 
surement®® induces a sizable error on the corresponding V -\- 6-jet cross sections (see 
e.g. Table 2). 

Beyond the parton level, the best theoretical predictions can be obtained via 
NLO QCD event generators (e.g. MCQNLD, POWHEG, SHERPA) that interface NLO 
QCD parton-level calculations with parton-shower Monte-Carlo event generators. 
In this context, care must be taken to assess and quantify the systematic uncer¬ 
tainties introduced by the parton-shower algorithms and by the hadronization and 
underlying-event models used. This is a non trivial task and need to be addressed on 
a process by process basis, taking into account the specifics of the generator used 
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and of the experimental set up. Nevertheless, general tests aimed at estimating 
the dependence of the theoretical prediction on parameters intrinsic to the shower¬ 
ing algorithms (like scale parameters other than renormalization and factorization 
scales) or the hadronization model should be always performed and the correspond¬ 
ing variations should be understood, confirmed by different approaches (i.e. codes), 
and included in the theoretical systematics. Several of the results that are collected 
in Tables 2 and 4 include an estimate of uncertainties from hadronization glob¬ 
ally quoted, together with an estimate of pile-up and underlying-event effects, as 
uncertainty from non-perturbative sources. 

In this context, it should also be mentioned that implementing 6 -quark initiated 
processes in NLO QCD Monte-Carlo event generators still presents some issues. We 
discussed in Section 2.1 how radiative corrections to tree-level 6 -initiated processes 
includes processes involving only light quarks in the initial state and b quarks in 
the final state (e.g.: gg —)■ bh^ is an 0{as) real-emission correction to bg —>■ 67 ). In 
order to automatically generate the radiative-correction processes from the corre¬ 
sponding tree-level process, the b quark should be consistently treated as massless 
or massive in both the initial and final state. Therefore, the traditional 5F prescrip¬ 
tion of considering a b quark as massless in 5-initiated processes and massive when 
it only enters as a final-state parton (the so called S-ACOT scheme^®) might be 
impractical, if not inconsistent, in the context of a NLO QCD event generator. Dif¬ 
ferent implementations of 5F in NLO parton-shower Monte Carlo have considered 
and partially adressed the problem with different methods, while a more rigorous 
solution is still being developed. 

To clarify this issue, let us observe that the treatment of processes initiated by 
heavy quarks and the rigorous definition of a heavy-quark parton density has been 
studied®^ and there is no conceptual obstacle in considering an initial-state heavy 
quark as massive. Indeed, a well-defined scheme to build fully-massive heavy-quark 
PDF have been proposed, the ACOT scheme,®^ and implemented in a few sets of 
PDF, namely CTEQ4-6(HQ). However, if the heavy-quark PDF is entirely gener¬ 
ated via perturbative evolution of gluon and light-quark density, it can be shown®® 
that the problem is equivalent to one in which the initial-state heavy quark is treated 
as massless and its mass is retained in processes that involve heavy quarks in the 
final state only. This is clearly very convenient in implementing a 5F calculation, 
since it allows to treat a large fraction of the contributing subchannels and their 
loop corrections using a massless heavy quark, greatly simplifying the complexity 
of the calculation. The corresponding heavy-quark PDF is defined in the S-ACOT 
scheme,®® where the heavy-quark mass (mb in our case) is kept only as a collinear 
regulator, and the corresponding collinear logs are resummed using the PDF evo¬ 
lution as reviewed in Section 2.1. Slightly different flavors of the S-ACOT scheme 
have been implemented in the CTEQ sets of PDF from version 6.5 on.®® Other sets 
of PDF (MSTW, NNPDF) implement slightly different matching schemes (TR®®’®^ 
for MSTW®® and FONNL®® for NNPDF®^) to account for the transition from above 
to below the heavy-quark mass thresholds, but all rely on the calculation of the b 
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initiated processes with rrib = 0 . 

Since a consistent implementation of ^-initiated processes in a parton-shower 
NLO Monte Carlo generator naturally call for treating the b quark as massive, 
one should probably consider adopting the original ACOT scheme and keeping 
the &-quark mass consistently throughout the calculation. Provided the radiation 
from initial-state massive quarks is correctly modeled in current NLO Monte-Carlo 
generators, this will allow accounting for all contributing sub-channels in each V + b- 
jet process, leaving the mutual balance between them to naturally adapt to the 
different energy regimes and kinematic regions. 

Let us conclude by adding that the incremental reduction of the current main 
sources of theoretical uncertainty will allow us to consider new level of precision 
in the future. For instance, there is no estimate at the moment of the effects due 
to higher-order electroweak (EW) corrections. They are expected to affect the to¬ 
tal cross sections by only a few percent, but they might have localized effects on 
distributions, and will have to be accounted for in the future. 


2.3. Photon isolation 

We would like to add in this review a special note on the issue of photon isolation, 
which concerns only the 7 -I- 16-jet and 7 -I- 26-jet processes. Since this issue is often 
discussed in the corresponding experimental analyses, it is important to clarify the 
difference yet compatibility between different possible theoretical approaches. We 
remind the reader that this has been thoroughly discussed in Ref. 9. 

Photons in a hadronic environment are usually distinguished into prompt pho¬ 
tons, when they are directly produced in the hard interaction, and secondary pho¬ 
tons, when they originate from the hadronization phase of a hadronic jet or the 
decay of unstable hadrons (e.g. 7 r° —)■ 77 ). While the production of prompt pho¬ 
tons can be described in perturbation theory, the production of secondary photons 
has associated non-perturbative effects that can only be modeled and can there¬ 
fore introduce a large parametric uncertainty in any given calculation. Since sec¬ 
ondary photons tend to preeminently occur in regions of the detector with abundant 
hadronic activity, in particular within or close to jets, their effect can be eliminated 
by imposing so-called isolation cuts which specihcally limit the hadronic activity 
around a given photon. Prompt photons become then isolated photons and can be 
easily disentangled. 

The main theoretical caveat in implementing a given prescription to isolate 
prompt photons from the hard interaction is that such procedure can veto regions of 
phase space responsible for soft QCD radiation and could therefore spoil the cancel¬ 
lation of infrared divergences between virtual and real corrections in a perturbative 
QCD calculation. As soon as some residual hadronic activity is admitted in the 
region around the photon, very energetic collinear final-state partons can produce 
a small parton-photon invariant mass and the corresponding collinear divergences 
need therefore to be reabsorbed into suitable non-perturbative fragmentation func- 
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tions (FF). The cross section for prompt-photon production is then given by, 

= cr2irect(M«)MF) + / dz^ ^iF, Mp), Mp), (5) 

Jo 


where represents the cross section for the direct-photon component while ai 

denotes the cross section for the production of a parton i that further fragments 
into a photon. The fragmentation of a parton i into a photon is represented by 
the corresponding photon FF, Di^^{z,Mp), where z is the fraction of the parton 
momentum that is carried by the photon, and Mp is the fragmentation scale. Ex¬ 
amples of available FF in the literature are by Bourhis, Fontannaz and Guillet (set 
I and II),®® and by Gehrmann-de Bidder and Glover.®® Fragmentation functions 
for final-state partons, like PDF for initial-state partons, are intrinsically non per¬ 
turbative and introduce into the calculation the same kind of uncertainty in the 
modeling of secondary photons that one originally wanted to eliminate. How rele¬ 
vant the contribution of fragmentation functions is depends on the chosen isolation 
prescription. 

Theoretical calculations normally use two main prescriptions denoted as fixed- 
cone and smooth-cone (or Frixione) prescriptions respectively. The fixed-cone pre¬ 
scription is commonly used in experiments and limits the hadronic activity inside a 
cone of radius Rq around the photon by imposing that the hadronic transverse en¬ 
ergy inside the cone does not exceed a maximum value, set by the experiment, 

i.e. 


^F;T(had) < , (6) 

eflo 

where Rq = , and A?; and Afi are the pseudorapidity and azimuthal- 

angle differences between the photon and any hadronic activity. After the isolation 
cut, the value of z is typically large, and since the FF are dominant in the low z 
region, the isolation procedure suppresses the fragmentation contribution substan¬ 
tially. 

Alternatively, the smooth-cone isolation prescription introduced by Frixione,^® 
limits the hadronic activity around a photon by imposing a threshold on the trans¬ 
verse hadronic energy within a cone about the photon that varies with the radial 
distance from the photon, i.e 

Eir 0{R - R^.-y) < cEf ( ^ ~ \ for R < Rg , (7) 

\ 1- COS Flq J 
i ^ ' 

where the i summation runs over all final-state partons in the process and E^ is 
the transverse energy of the parton (photon). Rq is the size of the isolation cone, e 
is an isolation parameter of 0(1), the exponent is normally set to n = I or 2, and 


i?i,7 = y^(T7^^TT(ATT ■ 
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The 0-function ensures that the i summation only receives contributions from par- 
tons that lie inside the isolation cone. R = Ri^j if there is only one parton inside 
the isolation cone. The r.h.s of Eq. 7 vanishes as R —t 0, thus the collinear configu¬ 
rations are suppressed while soft radiation is allowed to be present arbitrarily close 
to the photon. Since the collinear configurations are completely removed, there is 
no fragmentation component in Eq. 5. 

The fragmentation contribution in the pp{pp) —>■ bbj calculation is included at 
0{aa^). Due to the photon-isolation requirement, a photon cannot fragment from 
the tagged b/b quark. In the 7-1-26 case, the photon can only fragment off a light 
parton j, i.e. ai in Eq. 5 is the cross section for the pp{pp) —>■ bbj process calculated at 
LO, (cTi = aLo(pp(pp) bbj)). We notice that o lo{pp(.Pp) bbj) is finite since we 
impose a cut on the photon transverse momentum. For the 16-tag case, in addition 
to the same contribution present in the 26-tag case, the photon can also fragment 
off an unidentihed 6/6 quark. The LO pp{pp) —>■ bbj cross section is divergent in 
this case since the light parton in the final state can be soft and/or collinear. To 
overcome this problem, one should start from the pp{pp) —>■ 66 cross section at NLO 
in QCD. We also notice that when the photon is fragmented off of a 6/6 quark, 
terms proportional to ln{Mp/m'^), arising from the collinear configuration of the 
b ^ bj splitting in the (direct) pp(pp) —> bbj process, have to be subtracted to 
avoid double counting since those terms have been included and resummed in the 
6 quark-to-photon fragmentation function via DGLAP evolution equations. 

Having implemented both methods in 7 J- 26 production, very little if no dif¬ 
ference was found in Ref. 9. Similar conclusions have been reached in other studies 
involving the associated production of a photon and several jets.^^’^^ This should 
be emphasized to the benefit of the experimental measurements, since the exper¬ 
imental set up (where energy depositions are intrinsically discretized) only allows 
for a fixed-cone isolation criterion (while the smooth-cone criterion is a continuous 
prescription, even if discrete versions can be explored). 


3. Experimental Measurements and Theory Predictions 

Experimental measurements of vector-boson production in association with 6 jets 
started'^ in 2008^® and has been followed by more recent studies by all the major 
high-energy-physics experimental collaborations: 

• W + b jets by CDF,33 ATLAS, D 0 , 3 L 35 and CMS;^^ 

• Z + b jets by CDF,36 ATLAS,'*^’^^ and DO;^^ 

• j + b jets by D 0 , 39 Ao and CDF.^s 

More activity is expected in the near future, since, for instance, all the studies 
published so far by the LHC collaborations have been based on ^/s = 7 TeV data 
samples only. With the larger datasets collected in Run I of the LHC at ^/s = 8 TeV 


^We also note the earlier measurement in Ref. 73 of the cr[pp —> Z -f b]et)/a{pp —> Z -I-jet) ratio. 
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and expected in Run II at or above ^/s = 13 TeV, new possibilities will open and 
new challenges will be met, such that more accurate measurements will become 
available. 

One central and common feature of all measurements mentioned is the need for 
6 -jet tagging algorithms. The main characteristic for tagging 6 jets is of course the 
presence of a secondary vertex in the related jet. In order to increase the purity of the 
samples, other requirements can be imposed on, for instance, the tagging of decay 
products of the associated meson (typically leptons), or the impact parameters, 
invariant mass, and number of tracks associated with the secondary vertex. All 
these variables are combined in multivariate analyses. In the end the experiments 
report 6-jet tagging efficiencies between 35% and 50%, considerably lower than 6- 
tagging efficiencies around 85% reported, for instance, in studies of tt production 
(see for example Ref. 74). Indeed one must observe that, with the tagging algorithms 
employed in V -|-6-jet analyses, light-jet mistag rates are actually very low, at the per- 
mille level for light jets and at the few-percent level for c jets. This is necessary when 
studying R-|-6 jet signals, since much lower jet multiplicities are explored compared 
to tt studies. Unlike ti production, where one can also exploit kinematical constraints 
to exclude signals from light QCD jets, V + 6’s signatures can be overwhelmed by 
backgrounds like U -I- 1 light-jet with the light jet misidentified as a 6 jet. The net 
effect is that a high-purity 6-tagging algorithm has a reduced efficiency. 

An important issue in 6-jet tagging is the possibility of constraining the number 
of B hadrons in the tagged jet. Many of the algorithms employed so far look for jets 
that contain a B hadron, in principle associated with a 6 quark from the original 
hard interaction. But it has been observed that care must be taken with jets that 
contain pairs of B hadrons which can be associated with a gluon splitting into an 
almost collinear 66 pair. These jets are called merged 6 jets and denoted as (66) 
jets. Contributions of (66) jets in studies of signatures with 6 jets can be seen as 
reducible backgrounds. From the theoretical side, it is known that gluon splitting 
into heavy quarks needs special attention and improved modeling in parton show¬ 
ers (see discussion in Section 2). Strategies to reduce these backgrounds are then 
advantageous. In particular, studies like the one presented by ATLAS in Ref. 56, 
where techniques were developed to tag (66) jets, are very important. ATLAS in¬ 
deed found that they conld develop algorithms to reject (66) jets at the 95% level, 
while retaining 6 jets with a 50% efficiency. We notice that this latter study shows 
that a measurement on merged 6 jets is possible, a result that should be of great 
interest to the theory community that wants to improve the description of g —> 66 
splitting in parton showers. 

In the following subsections we review the most recent experimental measure¬ 
ments of U -I- 6 jets {V = W^, and highlight the comparisons to theoretical 

predictions that have been reported in the corresponding experimental papers. In 
Section 3.1 we present results for W + 6-jet production, in Section 3.2 for Z/y* -|- 6- 
jet production, and in Section 3.3 for 7 -|- 6 -jet production. We do not present all the 
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details of the experimental analyses, but limit ourselves to show tables and plots 
that well capture and summarize the main features of these studies with particu¬ 
lar emphasis on those aspects that have a direct influence on the comparison with 
theoretical predictions. In particular. Tables 1, 3, and 5 present a synopsis of the 
experimental set up for each analysis. In the first column we give the experiment 
(first the Tevatron experiments, CDF and DO, followed by the LHC experiments, 
ATLAS and CMS), the center-of-mass energy, the data set used, and the paper 
from where the results have been taken. In the second column we list the signatures 
that have been selected by each experiment, and in this context we denote by 6 a 6 
jet and by j a light jet, while we denote inclusive measurements (usually denoted 
in the experimental papers by at least one b jet or at least one b and one light jet 
or similar) adding a -\-X to the main jet signature. Finally, in the third colnmn we 
report the detailed kinematical cuts and vetos employed at the particle level, which 
should directly correspond to what has been used in obtaining the corresponding 
theoretical predictions. Experimental results and theoretical predictions for total 
cross sections (and ratios of) are listed in Tables 2 and 4. We have extracted the 
theoretical results from the corresponding experimental papers, and therefore they 
are not often directly comparable among themselves since they might correspond 
to different set ups and different calculations. We will highlight these differences in 
the discussion and emphasize what we can learn from the comparison. We will also 
include several figures from different studies to illustrate the comparison between ex¬ 
perimental measurements and theoretical predictions of total and differential cross 
sections 

3.1. Measurements of W hadronic production in association with 
b jets 

All studies presented for W 6-jet hadronic production consider leptonic decays of 
the W boson, employing either the electron channel (IF —)■ evg) or the muon channel 
(IF —)■ or both, while the decay into a tau lepton is considered as a background. 

Consequently the signatures studied contain missing transverse energy, due to 
the transverse momentum of the escaping neutrino, This is measured in 
the detectors by adding (vectorially in the transverse plane) energy depositions 
around the detector. In theory calculations such is identified with p^, which in 
experimental measurements can not directly be measured. In what follows we choose 
to always write while quoting kinematical cuts for the different experimental 
studies, as missing transverse energy is the fundamental object that can be studied 
experimentally.® 

The main backgrounds that affect the IF(—)■ Zt';) + 6-jet signatures are: IF + c jets 
and IF+light jets where one (or more) jet(s) is mistagged as a 6-jet; ti and single-top 


®We do this in spite of the fact that all the experimental studies use in constraining different 
signatures, because we think it is more appropriate. In fact, it would be better to identify with 
the negative sum of the momenta of the observed jets (light or heavy) in the transverse plane. 
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production; and multijet production (with some of the light jet mimicking leptons). 
Also Drell-Yan processes (with extra light jets) and diboson production can have 
relevant contributions as backgrounds. Kinematical cuts and vetos are often applied 
to minimize background, for example in extra leptons or jets to avoid contamination 
from top-pair production. 


Table 1. W &-jet experimental signatures and kinematics from all major high-energy- 
physics experiments. Details are presented in the text. 


Experiment 

Signatures 

Kinematics 

CDF,33 1.96 TeV, 1.9 fb"! 

W(^ lui) + b 

W{^ lui) + b + j 
W{^ liyi) + b + b 

pij, > 20 GeV, \n‘\ < 1.1, 

> 25 GeV 

Jets: 

Cone-based algorithm, R = 0.4 
> 20 GeV, |7?J| < 2 

DO,34 1.96 TeV, 6.1 fb”! 

W{^ lui) + b + X 

Muon channel: 

pi^ > 20 GeV, \r}f^ \ < 1.7, 

Electron channel: 

ptp > 20 GeV, 

|r?®| < 1.1, or 1.5 < |r?=| < 2.5 
> 25 GeV 

Jets: 

Midpoint algorithm, R = 0.5 
p*J. > 20 GeV, |»?*’| < 1.1 

ATLAS,42 7 TeV, 4.6 fb”! 

W{^ lui) + b 

Wi^ lui) + b + j 
W(-^ lui) + b + b 

pij, > 25 GeV, Ip‘1 < 2.5, 

> 25 GeV, Vfjf > 60 GeV 

Jets: 

Anti-/cx algorithm, R = 0.4 
p{ > 25 GeV, |p4j < 2.1 

AR{l,j) > 0.5 

CMS,43 7 TeV, 5.0 fb"! 

W (—>■ b b 

> 25 GeV, Ip^I < 2 . 1 , 

> 25 GeV 

Jets: 

Anti-fcx algorithm, R = 0.5 

> 25 GeV, |p*’| < 2.4 

Jet Veto: 

p^ > 25 GeV, |??4j < 4.5 


3.1.1. Experimental setups and total cross sections 

In Tables 1 and 2 we present a detailed account of the most recent studies of each of 
the major high-energy-physics experimental collaborations for W + 6-jets hadronic 
production. 

Experimental details are shown in Table 1. In the signature column we show the 
lepton channels employed, either IF(—>■ et'e) or IF(—>■ for the electron or muon 

channel respectively, and when these are combined we write IF(—>■ Ivi). Notice that 
tau decays are actually considered as background and no hadronic decay mode is 
considered since they are affected by large backgrounds. We pay special attention to 
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the jet multiplicity and, using the notation explained at the beginning of this sec¬ 
tion, we distinguish between exclusive and inclusive (• • • -I- X) measurement and we 
specify how many b jets and light jets are part of the signature. This is very impor¬ 
tant in understanding the comparison with theory. The kinematical cuts shown in 
Table 1 define the phase space considered to produce theoretical predictions at the 
particle level. We note that experimental analyses use more sophisticated selections, 
tailored to the optimization of the experimental apparatus and the available data, 
and subsequently transform these, via an unfolding procedure, to a set of simpler 
cuts and vetos to be used in theory Monte Carlo generators. Care must be taken 
to reduce as much as possible the unfolding procedure in order to minimize the 
sensitivity to the modeling of data in the Monte Carlo employed by the experiment. 


Table 2. W -h b-jet measurements by the major high-energy-physics experiments and corre¬ 
sponding theoretical predictions. The CDF measurement corresponds to a jet cross section, 
while all others are event cross sections. Details are presented in the text. 


Setup 

Cross sections in pb 

CDF33 

Experiment 2.74 di0.27(stat.) di 0.42(syst.) 

Theory 1.22 di0.14(scale) 

DO^-i 

Experiment {Muon channel) 1.04 ±0.05(stat.) ± 0.12(syst.) 

TYieory {Muon channel) 1.34 =1= 0.06(PDF)^Q'Q5(m^)) 

Experiment {Electron channel) 1.00 ±0.04(stat.) ± 0.12(syst.) 

TYieory {Electron channel) 1.28 ^q' 33 (scale) di 0.06(PDF)^Q'Q^('m^)) 

Experiment (Comfemeti) 1.05 di0.03(stat.) di 0.12(syst.) 

Theory (Combined) 1.34 (scale) ± 0.06(PDF)j:°;“®(mi,) 

ATLAS42 

Experiment (1 jet) 5.0 di0.5(stat.) di 1.2(syst.) 

Theorv (1 ietl 3 01 ±0.07(stat.)l°;g(scale) ± 0.04(PDF) 

Theory (1 jet) 3.01 ±o.08(NP)l“;« (DPI) 

Experiment (2 jets) 2.2 di0.2(stat.) di 0.5(syst.) 

Theory (2 jets) 1.69 ±0.06(stat.)+-° (scale) ± 0.04(PDF) 

^ ’ ±0.08(NP)ljj;Jg(DPI) 

Experiment (1 or 2 jets) 7.1 it0.5(stat.) it 1.4(syst.) 

rr,, ■ N . ±0.09(stat.)+°-®°(scale) ± 0.06(PDF) 

Theory (1 or 2 lets) 4.70 ) , 

^ ^ ±0.16(NP)l“;g(DPI) 

CMS-^^ 

Experiment 0.53 di0.05(stat.) di 0.09(syst.) di 0.06(theo.) di 0.01(lumi.) 

Theory 0.55 ±0.03(scale&PDF) ± O.Ol(NP) ± 0.05(DPI) 


The results of the measurements introduced in Table 1 are shown in Table 2, 
where we also collect some of the theoretical predictions quoted by the experimental 
papers. CDF reports jet cross sections, while all other experiments give event cross 
sections. We present the full set of reported uncertainties. For experimental results 
these are of statistical and systematic nature. Also, CMS separate the uncertainty 
due to Monte Carlo modeling in the extraction of selection efficiencies (“theo.”) and 
the one due to luminosity determination. Experimental statistical nncertainties are 
relatively low, of the order of 10%, for all measurements. It is to be expected that 
statistical errors at the total-cross-section level will be relatively marginal for all 
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future LHC studies (either at ^/s = 8 TeV or at -y/s = 13 TeV). The challenge for 
the experiments will then be to reduce systematic uncertainties (which so far vary 
between 10% and 30%) through a better understanding of their detectors. In par¬ 
ticular, reducing the systematic uncertainty associated to 5-jet tagging procedures 
might be crucial. For ATLAS and CMS, the larger luminosity environment of RUN 
II, with the larger number of pile-up events, will pose serious difficulties. 

The theory predictions listed in Table 2 are all parton-level predictions that 
include NLO QCD corrections as implemented in MCFM following Refs. 5, 27, 7, 
28. In the ATLAS and CMS studies, they have been corrected by non-perturbative 
effects. The theoretical uncertainty is separated in: dependence on the (unphysi¬ 
cal) factorization and renormalization scales (“scale”), dependence on the choice 
of parton distribution functions (“PDF”), dependence on the value of mb, non- 
perturbative corrections (“NP”), and dependence on the modeling of contributions 
from double-parton interactions (“DPI”). 

All theory results show a relatively large unphysical scale sensitivity around 
15% to 30%. Although a full calculation of Next-to-Next-to-Leading-Order (NNLO) 
QCD corrections might be far in the future, it might be possible to compute the 
gauge-invariant pieces most sensitive to scale variations. It has been noticed that 
this sensitivity comes mainly from the tree-like NLO QCD real corrections asso¬ 
ciated to a quark-gluon initial state. A step towards this has been achieved by 
computing I-loop amplitudes for Whhj,^^ and even more in Ref. 25 where full NLO 
QCD corrections have been computed. We notice that CMS, which has presented a 
measurement for the exclusive IF -I- 5-1-5 process, has estimated the scale sensitivity 
following the prescription of Ref. 75. This is a way to give a more meaningful theo¬ 
retical scale uncertainty since, for exclusive processes, the dependence on the scale 
variation tends to be relatively small (given the reduced impact from tree-like real 
contributions^’®). PDF uncertainties are mild and of the order of 2% to 4%. Non- 
perturbative uncertainties are mainly associated to hadronization and underlying- 
event models. ATLAS calculates them using the implementation of VF55 production 
in POWHEG BOX, which however is a purely 4F calculation and might not give the cor¬ 
rect information for W + h-\-X (while it would for IF-|-5-|-5or IU-I-5-I-5-I-A). CMS, 
on the other hand, estimates hadronization corrections comparing with a tree-level 
simulation obtained using MadGraphJ-PYTHIA. Since CMS measured a W+h-\-h cross 
section, using the aforementioned implementation in POWHEG BOX or in any other 
NLO parton-shower Monte Carlo generator (as mentioned in the introduction, Wbh 
is also available e.g. in MadGraph5_aMC0NLO) would have provided a more adequate 
information. In future predictions we cannot but recommend the systematic use 
of NLO QCD parton-shower Monte Carlo generators to directly assess the impact 
of these non-perturbative effects. Finally it is important to mention the large the¬ 
ory uncertainty quoted by ATLAS and CMS as due to double-parton interactions 
(DPI). ATLAS estimates that double-parton interaction has a large impact, of the 
order of 25%, on the total cross section, while CMS accounts for an effect at the 
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15-10% level. Even more, ATLAS shows that, when quoted differentially in 5-jet pt, 
these contribution mostly fall in the lower bins. It is interesting to notice that, as 
we will see in Section 3.2, the contribution of DPI to the -|- 5-jet cross section 
is estimated by CMS to be below 5%,"*® also occurring mainly in the lower bins 
in pip. This difference might be due to the better kinematical resolution achieved 
via the decaying products of the Zj^* boson {Zj^* —)• ll as opposed to IF —> Ui) 
which allows a better modeling of DPI. Given the theoretical challenge involved in 
rigorously describing DPI in QCD, more dedicated (data driven) studies are needed, 
in order to be able to reduce the DPI impact on IF + 5-jet measurements. 

Finally in Figure 3 we present a plot with a full account of the ATLAS mea¬ 
surements, including, on top of the NLO QCD predictions of Ref. 28 (MCFM 
4FNS+5FNS) other Monte-Carlo-based theory predictions: POWHEG+Pythia, 
based on the 4F calculation of IF55 but matched to a parton shower at NLO in 
QCD,^^ and ALPGEN-|-Herwig, matched to parton shower but without full NLO 
QCD corrections. Also, results are reported in jet-multiplicity bins and offer the 
possibility to assess separately the comparison between theory and experiments in 
each case. 

Looking at the results in Table 2 and Figure 3, one can see that the agreement 
between experiments and theory is relatively good, although within still sizable 
uncertainties. We notice that the agreement is definitely better for signatures with 
at least two jets. For single-5-jet signatures, CDF does have a large excess with 
significance around 3a. On the other hand, DO does not see such discrepancy, but it 
is important to notice that the measurement by DO is fully inclusive and it is an event 
cross section. Also, given the difference in their (both IR unsafe) jet algorithms, the 
comparison is difficult. On their side, ATLAS measurements are higher than the 
corresponding theory predictions, but with little significance (around 1.5cr). It will 
be interesting to see if such tension persists when ATLAS publishes updates with 
larger data sets. The comparison between different theoretical predictions given 
in Figure 3, with and without parton-shower and non-perturbative effects, shows 
clearly that different theoretical estimates are in relatively good agreement, giving 
confidence in particular in the quoted NLO uncertainty bands. Finally, we note 
that the slight experimental excess is shown clearly to come from the 1-jet bin. This 
seems to be confirmed by the complementary CMS measurement which focus on the 
IF+ 5 +5 signature, and find excellent agreement between theory and experiment. 


3.1.2. Differential cross sections 

In order to show how well experiments and theory compare over phase space, we 
show in Figure 4 a couple of plots from ATLAS^^ and DO.®® Notice that the differ¬ 
ential cross section from DO comes from a more recent stndy as compared to the 
one quoted in Table 1. 

Looking at the ATLAS differential cross section as a function of 5-jet pt it 
seems that the theory predictions deviate considerably from data at large transverse 
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Fig. 3. Total cross sections for W + fe-jet hadronic production measured by ATLAS"^^ and com¬ 
pared to theoretical predictions. 
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Fig. 4. Differential cross sections for various W + b-jet signatures, as a function of the b-jet px, 
measured by ATLAS^^ (l.h.s.) and DO®® (r.h.s.), and compared to theoretical predictions. 


momentum, where, however, uncertainties are also bigger. On the other hand, the 
DO spectrum seems systematically higher with respect to the theoretical predictions, 
although the ratio between theory and data seems constant. It will be interesting 
to study these differential distributions more in detail with future larger data sets, 
and assess the impact of DPI in various regions, as first studied by ATLAS,as 
well as the possibility that large bins might be “polluted” by (bb) jets, i.e. jets 
that contain a pair of B mesons. 
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3.2. Measurements of Zj')* hadronic production in association 
with b jets 

Measurements presented for Zj^* + 6 -jet production have employed both elec¬ 
tron and muon channels for the decay of the intermediate neutral gauge boson 
(Z/ 7 * —)• e'^ e~, p,'^ p,~). Therefore the main signal consists of a pair of opposite-sign 
leptons and associated b and light jets/ The main backgrounds in these studies are 
represented by top-quark pair production, W production with jets, Z production 
with light jets, inclusive Z production with Z decaying into tau leptons, diboson 
and single-top production. In order to enhance the signal, experimentalists choose 
kinematical cuts and apply vetos as needed. 


Table 3. + b-jet experimental signatures and kinematics from all major high-energy-physics 

experiments. Details are presented in the text. 


Experiment 

Signatures 

Kinematics 

CDF,36 1.96 TeV, 2 fb"! 

ir) + b + x 

76 GeV < Mil < 106 GeV 

Jets: 

Cone-based algorithm, /? = 0.7 
> 20 GeV, < 1.5 

DO,3'^ 1.96 TeV, 9.7 fb”! 

[z/T*(^U)+j+xJ 

70 GeV < Mu < 110 GeV 

Muon channel: 

'" > 15 GeV, '" > 10 GeV, 

IpA*! <2 

Electron channel: 
p?j. > 15 GeV, 

|??=| < 1.1, or 1.5 < |p®| < 2.5, 

Jets: 

Midpoint algorithm, R = 0.5 
> 20 GeV, lyi < 2.5 

ATLAS ,"‘6 7 TeV, 4.6 fb"! 

Z/y*(^ n) + b + x 

Zh* {-^ ll) + b + b + x 

plp > 20 GeV, Ip'I < 2.5, 

76 GeV < Mu < 106 GeV 

Jets: 

Anti-fcx algorithm, R = 0.4 
p{ > 20 GeV, |pl| < 2.4 

AR(l,j) > 0.5 

Cms, 48 7 TeV, 5.0 fb"! 

Z/7*(^ ll) + b 

Z/-i*{-^ ll) + b + X 
Z/7*(^. ll) + b + b + X 

plp > 20 GeV, Ip'I < 2.4, 

76 GeV < Mu < 106 GeV 

Jets: 

Anti-fcx algorithm, R = 0.5 
p{ > 25 GeV, < 2.1 

AR(l,j) > 0.5 


^This signature can be generated by both Z and 7 * although the first component is largely domi¬ 
nant. 
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3.2.1. Experimental setups and total cross sections 

In Table 3 we present details of the experimental setups employed on Z + 6-jet 
production studies, using the same format of Table 1. Unlike the case of IT -I- 6 - 
jet measurements, most of the results shown here are inclusive in the number of 
jets. Also, as it is needed in order to avoid the massless photon pole, a cut on the 
invariant mass of the lepton pair is always imposed around the Z peak. 


Table 4. Z/ 7 * + b-jet measurements by the major high-energy-physics experiments and corresponding 

theoretical predictions. CDF and DO measurements are ratio of cross sections. Details are presented in the 
text. 


Setup 

Cross sections or ratios 

CDF36 

Experiment (ratio to inclusive Z) (3.32 di0.53(stat.) di 0.42(syst.)) x 10“^^ 

Theory (ratio to inclusive Z) (2.3 — 2.8)(scale) x 10“^ 

DO^'^ 

Experiment (ratio) 0.0196 di0.0012(stat.) di 0.0013(syst.) 

Theory (ratio) 0.0206 ^Q'QQ^ 3 (scale&NP) 

ATLAS'^® 

Experiment [fb] (at least 1 jet) 4820 di60(stat.) ^ 3 gQ(syst.) 

Theory, MCFM [fb] (at least 1 fe jet) 5230 ±30(steF^ l^®°(scale&NP) 

Theory, 4FNS aMC@NLO [fb] (at least 1 jet) 3390 di20(stat.) ^^gQ(scale) 

Theory, 5FNS aMC®NLO [fb] (at least 1 jet) 4680 di40(stat.) ^ 3 gQ(scale) 

Experiment [fb] (at least 2 6 jets) 520 di20(stat.) 

Theory, MCFM [fb] (at least 2 6 jets) 410 ±10(sU^.^ lgo(scale&NP) 

Theory, 4FNS aMC@NLO [fb] (at least 2 6 jets) 485 di7(stat.) ^®Q(scale) 

Theory, 5FNS aMC@NLO [fb] (at least 2 6 jets) 314 di9(stat.) ^|Q(scale) 


Experiment [pb] (1 5 jet) 3.52 di0.02(stat.) di 0.20(syst.) 

Theory, MCFM [pb] (1 fe jet) 3.03 toJelscale) 

Theory, 4FNS aMC@NLO [pb] (1 b jet) 2.36 tolrlscale) 

Theory, 5FNS aMC@NLO [pb] (1 b jet) 3.70 to; 26 (scale) 

Experiment [pb] (at least 1 b jet) 3.88 di0.02(stat.) di 0.22(syst.) 

Theory, MCFM [pb] (at least 1 5 jet) 3.23 ^Q'^Q(scale) 

Theory, 4FNS aMC@NLO [pb] (at least 1 5 jet) 2.71 (scale) 

Theory, 5FNS aMC®NLO [pb] (at least 1 5 jet) 3.99 ^Q' 29 (scale) 

Experiment [pb] (at least 2 6 jets) 0.36 di0.01(stat.) di 0.07(syst.) 

Theory, MCFM [pb] (at least 2 5 jets) 0.29 ^Q'Q^(scale) 

Theory, 4FNS aMC@NLO [pb] (at least 2 5 jets) 0.35 ^Q'Qg(scale) 

Theory, 5FNS aMC@NLO [pb] (at least 2 5 jets) 0.29 ^Q'Q^(scale) 


In Table 4 we show the results of the experimental measurements, together with 
the theoretical predictions quoted by the experimental studies. Notice that, although 
CDF^® showed results for jet-level cross sections, we show only results for the ratio 
with respect to inclusive Z-boson production, since this is the observable that they 
compare to theoretical predictions. Also, as shown in Table 3, DO measurement is 
a ratio of event cross sections for inclusive production oi Z + b + X and Z + j + X. 
ATLAS and CMS results are event cross sections. For all measurements we present 
theoretical results based on parton-level calculations that include NLO QCD results. 
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Also, for ATLAS and CMS we include fully showered results including NLO QCD 
corrections for the hard interaction, both in the 4F and in the 5F schemes. 

Reported errors are of statistical and systematic nature for experimental mea¬ 
surements. For theory, reported errors correspond to statistical integration er¬ 
rors (“stat.”), sensitivity to unphysical renormalization and factorization scales 
(“scale”), non-perturbative corrections (“NP”), and choice of parton distribution 
functions (“PDF”). The latter has been estimated only for the ATLAS results, as 
they report theoretical predictions using three different sets of PDF (MSTW2008,^® 
CTIO,^^ and NNPDF23^®). 

The statistical errors of the experimental measurements are all below 6% except 
for the CDF result, which happens to use a fraction of their final data set. All 
systematic errors are relatively small, below 14% (except for the inclusive two b- 
jet measurement from CMS, which has 19%). Theoretical uncertainties tend to 
be dominated by a scale sensitivity of the order of 15%, with PDF uncertainties 
relatively low at the 5% level. ATLAS has included non-perturbative corrections to 
the parton level NLO QCD results due to DPI, and estimated them to be of order 
1%, much smaller than the corresponding corrections in W -1-6 jet production. As we 
commented in Section 3.1, this may be due to the better resolved kinematics that 
allows a better and more constrained modeling of Zj^* -\- 6-jets DPI contributions. 

Both ATLAS and CMS have included comparisons to theoretical results based 
on parton showers that consistently include NLO QCD corrections for the hard in¬ 
teraction. A summary at a glance of the ATLAS stndy for both Z/y* -|- 6 -|- A and 
Z/y* -1-26-1-A signatures is nicely illustrated in Figure 5. This is an important step, 
as we hope in the future to be able to consistently study all observables including 
quantum corrections fully showered and hadronized to reach hadron level predic¬ 
tions. Even more, when analyzing signatures with multiple jet bins, techniques for 
merging NLO QCD matrix elements for different multiplicities can be employed.^® 
Doing this one has a proper way of comparing theoretical results (with well esti¬ 
mated uncertainties) directly with the observables measured by the experiments. 
In particular, both ATLAS and CMS have employed the aMC@NLO framework. 
In the future, it will be relevant to systematically compare predictions using other 
NLO QCD parton-shower Monte Carlo, such as the POWHEG BOX and SHERPA. 

Comparing experimental measurements to the corresponding theoretical results 
obtained at NLO QCD or adding parton shower, in the 4F and 5F schemes, one 
can notice that shower and hadronization effects (not included in the CMS results 
from MCFM) have an impact of the order of 20%. Also, it is clear that the choice 
of scheme (either 4F or 5F) has a considerable impact at the parton-shower NLO 
QCD level. For example, it is clear that 5F results compare better with single-6-jet 
measurements, while 4F results compare best to two-6-jet measurements. 

Overall, at the level of total rates the Z/y*-|-6-jet signatures have good agreement 
between theory and experiments, at least when comparing to the proper ordering of 
the perturbative expansion (4F vs 5F). New studies with larger experimental data 
sets will help further testing this agreement, with particular attention to certain 
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features that appear over phase space as we will show in the following subsection. 
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Fig. 5. Total cross sections for hadronic production of Z/ 7 *+ > lb jet and Z/ 7 *+ > 2b jets from 
ATLAS^® compared to various theoretical predictions. See text for more details. 


3.2.2. Differential cross sections 

In Figure 6 we have included several plots of differential cross-section measurements. 
On the top left, we show a plot from CDF^® for Zjff + b-jet+X production as a 
function of the transverse momentum of the lepton pair (labeled as Px). On the 
top right we show a plot from ATLAS"^^ for Z/j* + 2 6-jet-I-A production as a 
function of the R separation between the b jets. Finally, on the bottom, we show 
two plots from CMS as published in Ref. 47. As CMS"^® did not show differential 
cross section, we have chosen to show the ones in CMS,"^^ which is a study of the 
associated production oi a Z boson and two B mesons. We do not include plots 
from because they only show ratios to Z/ 7 *-|-light jet production. 

In general the agreement between theory and experiments on all plots in Figure 6 
is relatively good. In the ATLAS plot for Zjff -|- 26 -|- A we notice the poor agree¬ 
ment between the 5F predictions and data, confirming what previously observed 
for the total cross section, i.e. that 5F is not expected to provide reliable results for 
processes with two or more 6 jets (notice that from a parton-shower point of view, 
this observable is computed at LO). Two of the plots in Figure 6 show distributions 
for Ai? (for 6-jets in the case of ATLAS, and B mesons in the case of CMS) and 
they both indicate that theoretical descriptions fail to describe the lower bins. A 
failure to describe the small-angle limit for the production of 6 jets (hadrons) could 
be ascribed to badly described gluon splitting into a slightly collinear 6-quark pair. 
One could then correlate these discrepancies with more collinear 66 splitting 
(producing (66) jets) which should show up in the tail of distributions in (and 
also Px^®*). Although statistics is poor, there is no evidence of this in the bottom 
right plot in Figure 6. Future studies might give more hints of correlations between 
discrepancies for low-angle distributions and large-px observables. In the end, this 










































June 25, 2015 1:20 WSPC/INSTRUCTION FILE Vbb'mpla 


Vector-boson production in association with b jets at Hadron Colliders 27 


w 0-3 


5 12 

I 

“o 8 


CDF Data 
■ MCFM Q^=m| 

. MCFM Q^=<p^ 


10 20 30 40 50 60 70 80 90 100 

[GeV/c] 


I Data 

l^i = 7 TeV, 4.6 ft’ - 
H MCFM 

‘^aMC@NL0 5FNS . 

aMC@NLO 4FNS - 
-T- ALPGEN+HJ 
-*-SHERPA 



1.2 

1 

0.8 




.1. 


.: 


. 

LL j 


. 

L 

: 

0.6 

0.4 





AR{b,b) 


- CMS, L = 5.2fb’ 

; = 7 TeV 

• data I 

■ MG 5F CTEQ6L1 I 

ut ] (pb) 

CMS, L=5.2 fb ' • data 

,(i = 7TeV ■ MG5FCTEQ6L1 

r allp^ 

♦ MG 4F MSTW2008 ~ 

A ALPGEN 4F CTEQ5M I 
T aMC@NLO MSTW2008 ^ 


♦ MG 4F MSTW2008 

mpilj A ALPGEM4FCTEQ5M 

J T T T T aMC(§>NLO MSTW2008 


(flr¥ ^ 

IQ-' 

Tjii : 

- f t 'I 

'. 

.. 1 1 ..~ 




§ 

1.5 

mi mi ittf ittj 


3.5 4 

ARo, 


p? , cut (GeV) 


Fig. 6 . Differential cross sections from CDF^® (as a function of the transverse momentum of the 
Z/ 7 * boson), ATLAS'^^ (as a function of the R separation between the two b jets) and 
(as a function of the R separation of the B hadrons and of the transverse momentum of the Zj'^* 
boson) for various Z/ 7 * + 6 -jet signatures compared to theoretical predictions. 


might point to a poor control of gluon splitting to massive quarks and require better 
theoretical modeling of (bb) jets. 


3.3. Measurements of 7 hadronic production in association with b 
jets 

Finally in this subsection we present measurements of the production of a photon 
7 in association with b jets. Unlike the analyses of the previous section, where the 
gauge bosons are tagged through their decay to leptons, the photon signatures are 
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recorded directly through activity in electromagnetic calorimeters. Pollution from 
signatures that mimic photons, like for example neutral pions decaying into pairs 
of photons in the electromagnetic calorimeter, have to be dealt with. We show next 
the measurements presented in the literature by CDF^® and D 0 .^° 


Table 5. 'y + b jets: signatures, kinematics. 


Experiment 

Signatures 

Kinematics 

CDF, 3 * 1.96 TeV, 9.1 fb"! 

'I + b + X 

Ef > 30 GeV, \y^\ < 1.04 

7 cone isolation {R = 0.4): 

< 2.0 GeV 

Jets: 

JETCLU alg., R = 0A 
plj, > 20 GeV, \y^ < 1.5 

Rj^j > 0.4 

00,-*° 1.96 TeV, 8.7 fb”! 

'I + b + X 
+ h + b + X 

pf > 30 GeV, lyTl < 1 
or 1.5 < \y~<\ < 2.5 

7 cone isolation {R = 0.4): 

Ei^° < 2.5 GeV 

Jets: 

Midpoint (Run II) alg., R = 0.5 
p{ > 15 GeV, \yf < 1.5 


3.3.1. Experimental setups and differential cross sections 

In Table 5 we present the experimental setups of the 7 + 6 -jet measurements. No¬ 
tice that both CDF and DO employ a cone isolation criterion to characterize the 
(prompt) photon. CDF has considered only the case of central photons (|j/'’'| < 1.4), 
while DO has presented results including both central photons (|j/'*'| < 1) and for¬ 
ward photons (1.5 < \y'^\ < 2.5). 

Both CDF and DO have shown their results differentially as a distribution in 
the transverse momentum of the photon (p^). We cannot therefore give a table of 
results, but will illustrate the comparison between experimental measurements and 
theory through the p^ distributions themselves. The comparison with the 'y + b + X 
measurement of CDF^® is illustrated in Figure 7 while the comparison with both 
the 'y + b + X and the j + b + b + X measurements of D0^®’‘^° is illustrated in 
Figures 8 and 9. For the comparison with the j + b + X measurements of both 
CDF^® and DO^® we show plots from Ref. 9 since they compare to the most up 
to date theoretical calculations, while for the comparison with the ^ + b + b + X 
measurements of DO^*^ we present plots taken from the experimental paper itself 
where results from Ref. 9 have been used. In Figure 9 data are also compared 
to results from a LO calculation with p^-dependent PDF,®°’®^ where however the 
uncertainty of the prediction is difficult to estimate and to compare with a full NLO 
calculation using integrated PDF. 

It can be appreciated that theory prediction in general describe the data. This is 
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a case even for 7 -I- 16 jet, where the 4F prediction is in better agreement with data 
than the 5F prediction since, as explained in Section 2, the dominant production 
mode at the Tevatron is via the qq channel, which is not influenced by large initial- 
state logarithms. Still, theory distributions seem slightly softer and fall a little below 
data at high p^. Measurements at the LHC will help pushing distributions to much 
higher and confirm or not these trends. In general there is considerable overlap 
between 4F and 5F predictions (denoted in the plots by 4FNS and 5FNS) for 1 6 -jet 
observables. As one might expect, from the missing extra 6 quark, 5F predictions 
tend to drop faster as a function of the photon transverse momentum. 
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4. Conclusions 

In this review we have presented the current status of electroweak vector-boson 
associated hadronic production with b jets. Measurements and theoretical predic¬ 
tions for IF-I- 6 -jet, Zjy -I- 6 -jet, and 7 -I- 6 -jet signatures were thoroughly reviewed. 
Results from the major high-energy-physics experimental collaborations, CDF and 
DO at the Tevatron, and ATLAS and CMS at the LHC, were summarized. 

By now we have reached a mature level in the study of these processes. This 
is very important given the impact that related studies can have in constraining 
couplings of the Higgs boson to 6 quarks, direct studies of heavy-quark parton 
distribution functions, and searches for new physics. 

Although in general we find good agreement between theoretical predictions 
and experimental measurements, some tensions are still present and important chal¬ 
lenges are left for the near future. With the larger data sets collected at ^/s = 8 TeV 
at the LHC, and with the even larger ones expected at and above y/s = 13 TeV, we 
expect improved statistical errors as well as reduced systematic errors (thanks to a 
better understanding of the detectors). Matching these reduced errors will be chal¬ 
lenging for theory predictions. Full parton-showered results, including NLO-QCD 
corrections and merging prescriptions for different jet multiplicities, as well as sys¬ 
tematic studies of hadronization and underlying-event models should become the 
standard approach. The need for NNLO QCD corrections as well as the inclusion of 
EW corrections and resummation of large logarithmic corrections could play a role 
in selected observables and will have to be investigated. More exclusive studies that 
include larger amount of heavy jets are now also theoretically and experimentally 
approachable and could lead to better control of the V -I- 6 jet processes at hadron 
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colliders. 
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